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Abstract 

We derive a lower bound on the ground state energy of the Hubbard model for 
given value of the total spin. In combination with the upper bound derived pre- 
viously by Giuliani our result proves that in the low density limit, the leading 
order correction compared to the ground state energy of a non-interacting lattice 
Fermi gas is given by ^'KaQuQd, where Qu(d) denotes the density of the spin-up 
(down) particles, and a is the scattering length of the contact interaction poten- 
tial. This result extends previous work on the corresponding continuum model 
to the lattice case. 



1 Introduction 

In recent years, much effort has been made to rigorously analyze the properties 
of dilute quantum gases at low temperature and low density. For repulsive pair 
interaction potentials, the leading order correction compared to the case of ideal 
quantum gases of the ground state energy and free energy of continuous quantum 
gases in the thermodynamic limit have been investigated in O [I3l [Ml [lOl [U [U [TBI [E] • 
In particular, in it was proved that the ground state energy per unit volume 
of a low-density spin 1/2 Fermi gas with repulsive pair interaction is (in units where 
h = 2m = 1) given by 

I (6^')'^' {q'J' + eT) + ^^<^QuQd + o{q^) . (1.1) 

Here, ^^^(rf) denotes the density of the spin-up (down) particles, and a > is the 
scattering length of the interaction potential. The total density of the system equals 
Q = Qu + Qd-, and {qu — Qdj/C^Q) is the average spin per particle. 
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The goal of this paper is to extend the analysis in [8j from the continuum to the 
lattice case. We restrict our attention to the case of a simple cubic, three-dimensional 
lattice. Without loss of generality, we choose units such that the spacing between 
two neighboring lattice sites is one; i.e., the configuration space for one (spinless) 
particle is . 

For simplicity, we consider the case where the interaction potential between the 
particles has zero range, i.e., only particles on the same lattice site interact. This 
is the simplest version of the Hubbard model, which was originally introduced as a 
highly simplified model for fermions with repulsive (Coulomb) interaction. For a 
review of the history, rigorous results and open problems, we refer to |18t [?]. 

An upper bound to the ground state energy of the desired form was already 
derived in hence we concentrate here on the lower bound. Our main result is 
given in Theorem [T] below. In combination, the two bounds show that the ground 
state energy per unit volume of the Hubbard model at low density Q = Qu + Qd and 
at given spin polarization {qu — Qd)/{'^Q) is given by 

eo(e'n, Qd) + 8TTaQu0d + o{q^) , 

where, as before, a denotes the (appropriately defined) scattering length of the in- 
teraction potential, and eo(^u, Qd) is the ground state energy per unit volume of the 
ideal lattice Fermi gas. 

1.1 Model and Main Result 

We consider particles hopping on the simple cubic lattice Z^. For spinless fermions, 
the appropriate Hilbert space TC{N) is the subspace of totally antisymmetric func- 
tions in L^(Z^^), with norm 

We define n{Nu,Nd) C L'^{Z^^-+^^'') as 

n{Nu,Nd)=n{Nu)^n{Nd). (1.2) 

Elements of 7i{Nu, N^) are thus functions of Ny_+Nc[ variables that are antisymmetric 
both in the first A^^ and the last A^^ variables. In appropriate units, the Hubbard 
Hamiltonian for spin-up particles and Nd spin-down fermions in a box [0, L]'^ is 
given by 

Nd Nu Nd 

1=1 i=l i=l j=l 

where xi, . . . x^^ and yi, ■ ■ ■ i/Nd ^^e the coordinates of the spin-up and spin-down 
particles, respectively. The usual lattice Laplacian is denoted by A; it acts as 
(A/)(x) = \x-y\=iif (y) ~ fi^))- We consider Dirichlet boundary conditions on 
[0,L], i.e., we restrict H to functions that vanish outside the cube [l,L — 1]^. The 



2 



thermodynamic limit corresponds to taking L — > oo, oo and — > oo in such 

a way that 

Qu = lim jj and Qd = lim -p- . 

Note that necessarily < gu(d) ^ 1 because of the antisymmetry of the wavefunc- 
tions. 

We remark that instead of considering two species of spinless fermions with par- 
ticle numbers Ny^ and A''^, as we do here, one can equivalently consider just one 
species of Nu + fermions with spin 1/2, and restrict to the subspace of total spin 
S = {Nu — Nd)/2. We will use the former formulation for convenience. 

The coupling constant g is assumed to be nonnegative and is allowed to take the 
value +00. The scattering length a of the interaction potential qSq^x is given by (c.f. 
m Eq. (1.6)]) 



g f 1 dk 



87r(57 + 1) ' 4Eii(l-cosfcO (27r)3 ' 

where k = {k^, k'^, k'^) £ R^. We refer to Section \T2\ for details. 
Our main result of this paper is the following. 

THEOREM 1. Let EoiNu,Nd,L) be the ground state energy of H in 173)) . // 

L ^ oo, Nu ^ oo, Nd ^ oo with Qu = limL-,00 Nu/L^ and Qd = ii^L~*oo Nd/L^ , 
then ^ 

liminf —rEo{Nu, Nd, L) > eoiQu, Qd) + S-KaQuQdC^ - o{l)) , (1.5) 

with 0(1) < C{ag^^'^)^^^^ for some constant C > 0. Here e(){Qu, Qd) is the ground 
state energy per unit volume of the ideal lattice Fermi gas, and q = Qu + Qd- 

As already pointed out in the Introduction, an upper bound of the desired form 
(jl.Sp was proved by Giuliani in [3], extending the method used in [8j. This shows 
that Eq. (jl.5p actually holds as an equality. 

It is easy to see that 

eoiQu, Qa) = \ [Q^'f [qT + qT) + 0{Q'h (1.6) 

for small q. Hence, at fixed a, the expression for the ground state energy of the 
continuous Fermi gas (jl.ip and the one for the lattice Fermi gas coincide up to terms 
lower of order o{q^). Theorem [T] is slightly stronger, however, since the error term 
is o(a^^), not o{q^). Hence (jl.Sp should be viewed as a result for small oq^^^ , which 
could be achieved either by making q small or by making a small. We point out, 
however, that the case of small a is much simpler than the case of small q, since the 
correct result can be obtained via first order perturbation theory in this case, while 
this is not possible for fixed a and small q. For small a and fixed q, the Hartree-Fock 
approximation becomes exact, as was shown in [1]. 

We state and prove Theorem [1] for the simple cubic lattice and zero-range in- 
teraction, for simplicity, but the method can be used in more general cases. For 



(1.4) 
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instance, different lattice structures could be considered, or longer-ranged hopping. 
Interactions of longer range could also be included, as long as they are non-negative 
and have a finite scattering length; and particles with more than two spin states 
could be considered, as in [8J. In combination with the methods developed in [16j . 
our technique can also be applied at non-zero temperature to estimate the pressure 
or the free energy. 

The proof of Theorem [T] follows closely the corresponding continuum result in 
[8], with several important and non-trivial modifications, however. One of the main 
ingredients is the generalized Dyson Lemma, stated in Lemma [U below, which allows 
for the replacement of the hard interaction potential g5xy by a softer and longer 
ranged potential, at the expense of the high momentum part of the kinetic energy. 
The proof of the corresponding Lemma 4 in [8j uses rotational invariance of M'^ in 
an essential way, and hence does not extend to the lattice case of Z^, where such a 
symmetry is absent. (See also the discussion in [3]). Our new Lemma [1] does not 
rely on this symmetry, however. 

Another important estimate in [8j that does not carry over to the lattice case is 
a bound on the average number of particles that are close to their nearest neighbor. 
The estimate in [S, Lemma 6] uses an inequality by Lieb and Yau \12\ Thm. 5] whose 
proof also relies on the rotational invariance of M^. In Lemma[l]below we will present 
a weaker version of this inequality which is equally valid in the lattice case. 

The strategy of the proof of Theorem [1] is similar to the proof in [8]. First, we 
separate the Hamiltonian H into two parts, the low momentum part of the kinetic 
energy on the one hand, and the high momentum part of the kinetic energy together 
with the interaction energy on the other hand. The first part is larger than eo(^u, Qd), 
while the second part can be bounded from below by a softer interaction potential; 
this is the content of Lemma [TJ In Lemma [21 we shall show that at low density the 
one particle density matrix of the ground state of the Hubbard model is close to a 
projection, namely the projection onto the Fermi sea. This information will then be 
used to prove that the expectation value of the softer interaction potential is given 
by SiraQuQd to leading order. To bound some of the error terms, we need a bound 
on the expected number of particles whose distance to the nearest neighbor is small. 
This will be accomplished in Lemma [3l 

In the next section, we shall state some preliminaries and introduce the notation 
used throughout the proof. Section [3] contains the main three Lemmas, which we 
have already referred to above. Finally, in Section |4] the proof of Theorem [T] will be 
given. 

2 Preliminaries 
2.1 Notation 

We start by introducing some notation that will be used throughout the proof. First, 
the gradient operator V = (V"*^, V^, V^) on L^(Z^) is defined as usual as 

(V7)(x) = /(x + e*)-/(x), 
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with denoting the unit vector in the i'th coordinate direction. Its adjoint is given 
by (V't/)(a;) = f{x - e') - f{x). The Laplacian can then be expressed in terms of 
the gradient as 

-A = • V = V • vt . 

For any subset A C Z^, we denote by 9a its characteristic function. It will be 
convenient to introduce the notation 



Note that this is a nonnegative operator which plays the role of the (Neumann) 
Laplacian on A. For A C 7,^ bounded, we denote by Pa the projection onto the 
normalized constant function on A, 



\0a){0 



— un 11 2 



II^A||2 

For h G L^(Z^), we define the convolution operator as 

{CJ}){x) = h* tj;{x) =^h{x- y)tp{y) . 



Its adjoint is given by (QV')(a;) = J2yez^ - x)ilj{y). 

We recall also the natural definition of the Fourier transform, mapping L^(Z^) 
to L2([-7r,7r]3). For p= {p^,p'^,p^) G M■^ \p'\ < tt, 

Its inverse is given by 

{2tt)'^ 7[-,r,7r]3 

Using the above definitions, the following properties are easily verified: 

S^oip) = 1 (2.1) 
h * tp{p) = h{p)^(j)) 

i2 _ / i.-7r/„M2 dp 



„r„2- / mP)\ ...3 



In particular, if M and M' are two functions satisfying 

|M(p)P + |M'(p)|^ = 1 for all p e [-tt, 7r]=^ , 
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we can decompose the Laplacian A as 



(2.2) 



We will use this decomposition in the proof of Theorem [T] in order to separate the 
kinetic energy into the high momentum and the low momentum parts. 
Finally, it will be convenient to introduce the operator 



It has the property that 



(2.3) 



(2.4) 



xeA,y(fA 
\x-y\=l 



2.2 Scattering Length 

We denote by ip the solution of the zero-energy scattering equation 

- Aip{x) + ^g6o^xV>{^) = 
with boundary condition lim^x\^oo ^{x) = 1. It is given by [3l Eq. (1.5)] 

e^P-^ dp 



tf{x) 



1 - Ana 



[-n,7r]^ 2X:?=i(l-cospJ) (2vr)2 



(2.5) 



(2.6) 



where a is the scattering length (|1.4p . It can be shown [15] that there is a constant 
C > such that 



ip{x) - 1 + rr 



13 • 



(2.7) 



Note that, in particular, a = lim^x^^^{(p{x) — l)\x\. It can be readily checked 



that 



Aifix) = lg^iO)=4TTa. 



Another property of if we will need is [31 Eq. (1.7)] 



\x-y\=l 

for any simply connected domain A containing the origin. 



(2. 



(2. 
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2.3 Non-interacting Fermions 

We recall here briefly the ground state energy of non-interacting (spinless) fermions 
on the lattice 1? . For "p G [— 7r,7r]^, the dispersion relation will be denoted by 



E(j)) = 2^ (1 - cosp^) . 

i=\ 

For given density < ^ < 1, the Fermi energy Ef{Q) is determined by 

(27r)"3 / dp = g. (2.10) 



'E{p)<Eiig) 

The ground state energy per unit volume in the thermodynamic limit is then 

e(Q) = {27T)-^ [ E{p)dp. 
JE(p)<Ei{e) 

For spin 1/2 particles with spin-up density Qu and spin-down density Qd, the ground 
state energy is thus eo(£'n, Qd) = aiOu) + aiOd)- 

3 Auxiliary Lemmas 
3.1 Lemma One 

As mentioned in the Introduction, Lemma [1] is the main tools of this paper. It is 
similar to Lemma 4 in [8j . This lemma allows for bounding the hard interaction gSo^x 
from below by a softer interaction at the expense of the high momentum part of the 
kinetic energy and some error terms. 

LEMMA 1. For r E N, let A{r) denote the cube A{r) = [-r,rf n Z^. For any 
function h G L^(l/') satisfying 1 > h{p) > 0, let 

fr{x)= max \h! [y) — h! {x)\ , where = 1 — . (3-1) 

yex+A(r) 



U = i2R + l)-^eAiR), (3.2) 



For R gN, let U, W and V denote the nonnegative operators 

W = WirfR fR{x) (3.3) 

and 

V = {2R + l)-^[eAiR)-PAiR)]- (3.4) 
There exists a constant C > such that for any R> C and 0<e<l, 0<7?<1, 

w cv 



Cl[Vt.0^(^)V],C^ + |4,o>47ra 



{l-e){l-rj)U 

e rj 



(3.5) 
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Compared with the result in [HJ Lemma 4], there is an additional error term 
V on the right side of (j3.5p . There is no restriction that U has to vanish at the 
origin, however, which was necessary in [8]. We note that the norm of U is given by 
||?7|| = {2R + which is much smaller than the norm of {g/2)6xfl (which is g/2) 
for our choice of i? ^ 1 below. 

Proof. We are actually going to prove the stronger statement that 

(ct[Vt • eV].C7, + f 4,o) > 4™ ((1 - 5)(1 -rj)U---—) (3.6) 

for any tp G L^(Z^). Here and in the following, we use the shorthand notation 
( •)^ = ■ {ip). The non-negative function is defined by 

^ R~l 
r=R 

where we denote by R the largest integer less than R/2. Since < O^^ji-j, (13. 6p 
implies (|3.5p . 

To prove (13. 6p . we first define Br as 

5R = (v|cI[vt.ev],|<^) + (^||5.,o|<^). 

Here, (p is given in (|2.6p . Using Schwarz's inequality, if follows that |-B_rP is bounded 
from above as 

< (c^[vt . ev],c^ + |4,o)^ ([v^ • ev], + |5,,o)^ • (3.7) 

By the definition of ip, ([V^ • GV]s + f i^^.o) < ("■^ + 2^x,o) = ^vra. Hence we see 
that the left side of (13. 6p can be bounded from below as 

{cJ[Vt.eVl,C, + |*,„)_^>lM. ,3,8) 

Define x G L'^i^^) via 

Ix) = [Vt . eV],|<^) + A|v9) = [Vt . QVU^) + |4,o|</') . 

Alternatively, using (j2.3p . |x) = {R — R)^^ J2flfi'^A{r)\^)- Hence x is supported in 
A{R) \ A{R — 2). Moreover, x(^) is a non-negative function for R large enough, as 
can be seen from ()2.4p and the asymptotic behavior of ip in (j2.7p . 

Let also |a) = [V^ • 0V]s|99). Since a{x) = x{x) ~ (5/2)^3;, o<^(0), also a is 
supported on A{R). Moreover, ^^a(x) = because of (|2.8p and (j2.9p . Since 

i?ii = (V'|x)-(V'|Ci,|a). 
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Recall that a{x) = 0, and a is supported on A{R). Hence 



(C^,a)(x) 



^ (^h'{y - x) - h'{x)^ a{y) 



with //J defined in (j3.ip . Here we used the fact that is non-negative and sup- 
ported away from the origin, hence = J2x xi^) + 5'V(0)/2 = 9^(0) = Sna. 
In particular, we conclude that 

\Br\ > KV|x)|-8™ 
Using Schwarz's inequality and the definition of W in (|3.3p . we thus obtain 



\BR\'>{l-e){ij\x){x\i^) 



(3.9) 



To get a lower bound on we use again Schwarz's inequality, as well as the fact 

that X is supported on A{R). We obtain, for < r/ < 1, 

|X)(X| > (1 - v)PAiR)\x){x\PA(R) + (1 - V-') {Oa(R) - Pa(r))\x){x\{Oa(R) - Pa{R)) 
> (1 - ^){A7,afj0^ + WxWl (1 - 7?-^) {Oa^r) - PAiR)) ■ (3.10) 

Here, we have again used the fact that x{x) = 47ra. 

To conclude the proof, we have to show that Hxlli — const, a^/i?^ for large R. 
This follows from the fact that 



X{x) < 



sup \ip(x) - (p(x + e)\ < Ca/R^ 

ri- K \x\>R,\e\=l 



for large i?, as can be seen from the asymptotic behavior (j2.7p . Inserting the inequal- 
ities (I3.10p and ()3.9p into ()3.8p . we arrive at the desired result ()3.6p . ■ 

If 2/2! > 2\/3i?, the cubes of side length 2R centered at yi and 1/2, respectively, 
are disjoint. Hence we can obtain the following corollary of Lemma [TJ 

COROLLARY 1. Let U, W and V be as in LemmaUi and let Uy = TyUTl, Wy = 
TyUTy and Vy = TyVTy, where Ty is the translation operator {Tyil)){x) = '4j{x — y). 
Ifyi,---yn satisfy \yi -yj\ > 2^/3R for all i 7^ j, then 



i=l 



i=l 



{l-e){l-rj)Uy^ 



(3.11) 
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3.2 Lemma Two 

Recall that the Hilbert space n{Nu, N^) is the subspace of L^{I?^^"+^-i'^) of functions 
that are separately antisymmetric in the Nu spin-up variables and the spin-down 
variables. For <I> G Ti.{Nu, Nd), let 'ju and 7^ denote the reduced one-particle density 
matrices of ^ for the spin-up and spin-down particles, respectively, with Tr'ju = 
and Trjd = N^. 

For m G Z^, we define the functions fm{x) G L?'{'L'^) as 

fra{x) = {L + l)-3/2 exp [2^im ■ x{L + l)-^) e^oM^) ■ (3-12) 

Note that the /m's are orthonormal functions. For any function ip supported on 
[1, L — 1]^, the expectation value {tpl — A|^) can be expressed as 

(V'|-A|V>= Yl E{27rm/{L + 1))MU)\\ (3.13) 

me[-L/2,(L+l)/2]3 

For Men, let ^(M) denote the projection 

E \fm){fm\. (3.14) 
E{2iTm/{L+l))<E{{M/{L+iyi) 

It is easy to see that 

lim M-^TrC{M) = 1 (3.15) 

M— >oo,L^oo 

in the thermodynamic limit L ^ 00, M ^ 00 with M/L^ g for some < ^ < 1. 

Let C(iVu) = and ^{Nii) = Cd for simplicity. As mentioned in the Introduction, 
we shall show in Lemma [2] that the reduced one particle density matrix juid) of ^ 
ground state $ of is close to S,u{d) iii an appropriate sense. 

LEMMA 2. Let <I> G T~(-Nu,Na- Assume that, in the thermodynamic limit N^i^^-j — > 00, 
L ^ 00 with = iV„(rf)/L3 fixed, 



lim sup 



E^'iAx, - E^JAy. ^ ) < eoiQu, Qd) + Cag' (3.16) 



for some C > independent of a and g. Then 



limsup --^Tr[7„(rf)(l - iu(d))] < const, g^fag^. (3.17) 



Proof. The proof is parallel to the proof of Lemma 5 in [8j, following an argument 
in m. ■ 
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3.3 Lemma Three 

For given points yi ■ ■ ■ i/Nd ^ ^'^) let I_R(yi, UNa) denote the number of yi's whose 
distance to the nearest neighbor is less or equal to 2^/3R. Because Corollary [1] can 
be applied only to those T/j's that stay away a distance larger than 2\/^R from all the 
other particles, we will need an upper bound on the expectation value of Ir in the 
ground state of the Hubbard Hamiltonian. This will be accomplished in Lemma E] 
below. It states that as long as R is much less than the average particle distance 
Q~^^^, the expectation value of Ir in the ground state of H is small compared to the 
total number of particles N. 

LEMMA 3. Let $ G 7i{Nu, N^). Assume that for some constant C > 0, indepen- 
dent of Qu and Qd, 



N ' 



T.^JiA, <Cg'^'-'. (3.18) 



-,2/3 



Then 

{'^\lR{yi,...,yNdm < const. N{{R + if 0f/\ (3.19) 

Proof. With Di denoting the distance of yi to the nearest neighbor among the y^ 
with /c / i, we can write 

iR{yu ■ ■ ■,yNd) = ^^[o,2V3R](^i) ■ 

1=1 

Here, 0\n2^R] denotes the characteristic function of [0,2\/^R]. It follows from 
Lemma In below that 

Nd 211/2 
i=l i=l x€Z'i 

for any 6 > 0. Using the assumption ()3.18p and optimizing over the choice of b, we 
arrive at the result. ■ 

It remains to prove the bound (j3.20p . which is a special case of the following 
Lemma. Its proof uses a similar decomposition method as in |llj . 

LEMMA 4. Let f be nonnegative, with J2xez^ fi\^\) ^ '^^^ denote the 

distance of Xi to the nearest neighbor among all points Xj with j ^ i. On the sub space 
of antisymmetric N -particle wavefunctions in L'^{T?^), 

N 011/2 

(-A. - /(A)) > --^N f{\x\f^ . (3.21) 
j=i zez3 

Proof. Let A^i be the largest integer less or equal to N/2, and let = N — Ni. 
Consider a partition P = (vri,7r2) of the integers 1,...,A^ into two disjoint subsets 
with A^i integers in vri and integers in tt2. For a given P and « G vri, we define 

D[ = min{|xi -Xj\,j e 7T2} ■ 
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It is easy to see that 



i=l \ i/ p jg^^ 



where the sum runs over all (^_^) partitions of {1, ... , N}. This follows from the fact 
that for given i and j, the probability that a partition P has the property that i G tti 
and j G 112 equals NiN2/{N{N - 1)) > 1/4. 
In particular, we see that 



i=i ^ p ie-Ki ^ ^ 

For fixed Xj, j & 1^2, we can use the Lieb-Thirring estimate [SJ Thm. 5.3] to conclude 
that 

Note that the antisymmetry of the wavefunctions is essential here. We can estimate 

f{D^r^<N2Y.f{\x\f\ 

Using in addition that n'^^"^ N2 < , we arrive at the statement. ■ 

4 Proof of Theorem [1] 

We write the Hamiltonian H in (11.30 as 



I Nu Nu Na \ / Na Nu N^ \ 

H={-E^^^ + yEE ^^^^y. + - ^ -^^^ + ^5 E E • (4-1) 
\ i=i i=i j=i / \ i=i i=i j=i / 

Recall that we restrict H to functions that are antisymmetric both in the x and the 
y variables, and are supported on the cube [1,L — ij^C^u+^d) jj^ the following, we 
are going to derive a lower bound only on the first term. The lower bound on the 
second term can be obtained in the same way by simply exchanging the role of x and 
y. Our bound is not a bound on the ground state energy of this first term, however, 
but rather estimates the expectation value of this term in the ground state of the 
full Hamiltonian H. 

First, as mentioned in the Introduction, we decompose the kinetic energy —A 
into a high and a low momentum part. Let again E{{g) denote the Fermi energy of 
an ideal gas of spinless fermions at density g, defined in (j2.10p , and let 
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Here, [ • ]+ = max{ • , 0} denotes the positive part. Moreover, let M'{p) = y 1 — M{p)'^. 
As pointed out in Eq. ()2.2p . we can decompose the Laplacian as 



We first claim that 



lim — rinfspec 



1=1 



(4.2) 



The proof follows in exactly the same way as the proof of Eq. (64) in [8j, using an 
argument in [6]. 

We proceed with the high-momentum part. Let / : ^ be a smooth, radial, 
positive function, with l{p) = for \p\ < 1, l{p) = 1 for \p\ > 2, and < l{p) < 1 
in-between. As in |8J, we choose hs{p) as 



hsip) = Ksp) ■ 
Since hs{p) = for \p\ < 1/s, we can estimate 

M{pf 





> 


[ E{p) \ 


+ 



E{{Qu) 



mm 



bl>i/s 



E{p) 



(4.3) 



(4.4) 



(Here, the minimum is taken over p G [— vr, vr]'^, \p\ > 1/s.) In particular, this implies 
that 



1 



E{{0u) 



min|p|>i/, E{p) 

Since E{p) ~ |pp for small \p\, and Ef{gu) < const, gi^^ , we obtain 



EfiQu) 



mm\p\>i/sE{p) 



> 1 — const, s 



2 2/3 



(4.5) 



(4.6) 



as long as s » 1. 

We note that with this choice of hs[p) the corresponding h'g{p) = 1 — hs{p) is 
a smooth function that is supported in \p\ < 2s"^. As in [8j, we conclude that the 
corresponding potential W{x) defined in Lemma [T] satisfies (for \ < R< const, s) 



R^ 

W{x) < const. — and W{x) < const. — 



(4.7) 



for some constants depending only on the choice of I. Moreover, if \yi — yj\ > 2^/3R 
for all i ^ j, then 

Wnix - y.i) < const. (4.8) 

i=l 

independent of x and A'^^. 
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We will now use Corollary [T] to get a lower bound on the sum of the high mo- 
mentum part of the kinetic energy and the interaction energy. In order to be able 
to apply this corollary, we have to neglect the interaction of the x particles with 
those y particles that are not at least a distance 2\/3R from the other y particles. 
Let us denote Y = (7/1, . . . ,yNa)j and let Y' be the subset of Y containing those yi 
whose distance to all the other 7/j's is larger then 2\/3i2. Note that, by definition, 
the cardinality of Y' is \Y'\ = Nd — /^(yi, • • • ,yNd)j with defined in Section [331 
Moreover, let Y" C Y' be the set of yj E Y' whose distance to the boundary of 
[0, L]^ is at least R + l. As argued in [8], \Y"\ > \Y'\ — const. {L/Kf and hence, in 
particular, 

lim L-^{^\{\Y'\ - \Y"\)\ = (4.9) 

L— >oo 

in the ground state $ of H. 

Applying Corollary [H together with (j4.5p . we obtain that for a given configuration 

of y. 



i=l 

with Wy defined as 
Wy = 



Ei{Qu) 



min|p|>i/s^(p) 



Nu 
+ i=l 



(4.10) 



^ 4vra((l-e)(l-r?)C/, 



cv, 



Vj 



(4.11) 



For any $ G 7i{Nu, N^,), we can express the expectation value of X^Jify]?, as 



Nu 



[wy] 



i=l 



(4.12) 



Y 



where we denote by qy the distribution function oiY = [yi ■ ■ ■ yN^)-, that is, 

QY = m5Y={y,...y^^}\^) 

and 7y denotes the one-particle density matrix of $ for fixed Y , i.e., 



1y{x,x') 



QY 



^ ^{x,X2, . . . ,XNu,y)^{x' ,X2, . . .,XNu,yT 



X2,---,XNu 



Note that < 7y < 1 and Tr7y = N^- Moreover, qy = 1 and qy'Jy = 7m, 
which is the one-particle density matrix for the spin-up particles introduced earlier 
in Section [3121 

Prom now on, we will consider <^ to be a ground state of H. The assumptions 
of Lemma [2] and [3] are clearly satisfied for this as the upper bound derived in [3] 
shows. 

Recall that (,u denotes the projection £,i^u) defined in (j3.14p . We write wy = 
wy,+ — Wy,-, where wy,+ > stands for the part of wy in (j4.1ip containing U, 
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whereas wy- > is the part of wy containing W and V. As proved in [HI Sect. V.C] 
we have, for any 6 > 0, 

Tr[jywy] > TV[e„u;y,+](l - <5) - Tr[e„u;y,-](1 + S) (4.13) 
- (1 + 6-') i\\wy,+ \\ + llw,-!!) T¥[7y(l - Cu)] - \\wy\\Tr[Ul - 1y)] ■ 

We are now going to bound the various terms on the right side of (|4.13p . First, 
using X^^jg^s U{x) = 1 and the fact that has a constant density, we have 

Tr[CunJy,+] = J^^i^ " " r/)4vra|y"| . 

To estimate hmx^oo L""^ J2y QY'^r[CuWy^+], we can use (14. 9p . \Y'\ = N^ — Ir{Y) and 
Lemma [3] to conclude that 



hm L'^ 0yTT[^uWy,+] > (1 - e)(l - ri)A-KaQuQd - const. a£i^((i? + 1)^^)^/^ 



Here, we have also used that limj;,_^oo L ^Tr[^u] = g^, which follows from (j3.15p . 
Analogously, using (j4.7p and the fact that for any ijj and fixed y. 



lnlV^) = ^.^An6 E IV'(x)-V'(x')l'<^^ max 

z(^zJx-t- ij" — z (x,x')ej/+A(_R) 



(x,x')e?/+A{i?) 
it is easy to get the upper bound 

■i?2 2/3- 



(^2 
2 + 



?7 



Moreover, using (j4.8p and the fact that the distance between two yj £ Y" is at least 
2\/^R, as well as r/ < 1, we find that, 

llwll < + \\wy,-\\ < const, a (^^-J^ + . (4.14) 

The bound in Lemma [2] implies that 
hm V ^>yTr[7y(l - ^J] = lim L-3Tr[7,(l - ^u)] < const. Q{a^g)^/^ . 

L— >oo ^— ' L— >oo 

Finally, the last term in (j4.13p can be bounded as 

IToo^'^ E QY\\wy\\Tr[^u{l - 7y)] < const, a + £^(a3£,)i/6 , 

where we have used (|4.14p as well as the fact that Tr[^u(l — 7y)] = Tr[7y(l — ^„)] + 
TrKn ~ 7y]- The last term, when averaged over Y, is o{N) in the thermodynamic 
limit, i.e., J2y QyT^[Cu - 7y] = o{N). 
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Collecting all the bounds, and applying the same arguments also to the second 
term in (j4.ip . we arrive at the lower bound 



lim ^Eo{Nu,Nd,L) 
> eoieu, Qd) + STTaQuQd 



I 2 2/3 R^q"^^ 

l-e-ri-5 - const. ( s'^g'^''^ + — + 



const. [ {R^gf/^ + ( + 



3. 



f?)l/6 / 1 1 



TfR^g es'^Rg 

Here, we have assumed that i? > 1 and that s ^ 1 in order to be able to apply the 
estimate (14.61). If we choose 



R = ^,-l/3(„3^)l/30 ^ ^ ^ ^-1/3(^3^)1/90 ^ ^ ^ 5 ^ ^ ^ {a'g)^^^ (4.15) 

this implies that 

^irn^ ■^Eo{Nu,Nd,L) > eo{gu, Qd) + STraguQd (l - const, {ag^^^)^^^^^ , 

which is the result stated in (jl.Sp . 

Recall that in order to be able to apply Lemma [H it is necessary that R> C for 
some constant C > 0, which for our choice of R in (I4.15|) is the case if 

{ag'/^f>Cg'/'. (4.16) 

Under this assumption, also s ^ 1 is satisfied for small a'^g. For fixed a, (j4.16p holds 
for small enough g. However, if a is very small, (j4.16p is violated. In this case, one 
can obtain our main result (II. Sp actually much easier. One simply omits the use 
of Lemma [U altogether, and applies our perturbative estimate ()4.13p directly to the 
interaction potential {g/2)5xy Notice that for small a, g ^ Svra is also small. The 
resulting bound is 



lim -^EoiNu, Nd, L) > eo{gu, Qd) + gQuQd 



— const. 

og 



(4.17) 



Choosing 5 = g~^^'^{a? g)^/^"^ and noting that g > Sna and g > const, (a'^g)^^^^ in the 
parameter regime we are interested in here, this yields (jl.Sp in case (j4.16p is violated. 
This finishes the proof of Theorem [TJ 
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